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The Runll physics program at the Tevatron started in spring 2001 with protons and antipro- 
tons colliding at an energy of y / s=1.96 TeV, and it is carrying on with more than 500 pb _1 of 
data as collected by both the CDF and D0 experiments. Recent results on beauty production 
cross section measurements are here reported. 



1 Introduction 

Beauty production measurements provide precision tests of perturbative QCD. At the Tevatron, 
the total bb cross section is about 50 fib, resulting in an event rate of few kHz, thus, in a very high 
statistics. In the past years, the publication of " controversial" beauty cross section measurements 
by CDF and D0 in RunI, has led to many developments both in the theoretical calculations 
beyond NLO and the experimental approach, resulting in a better agreement between data and 
theory as it is shown in this paper. A major role, concerning theoretical improvements, has been 
played by the implementation of the so called Fixed-Order with Next-to-Leading-Log (FONLL) 
calculation^. Also, there have been substantial changes in bottom fragmantation function as 
extracted from experimental data and for the used Parton Distribution Functions (CTEQ6M). 
On the other hand, from the experimental side, there have been many improvements in the 
treatment of data inputs, for instance avoiding deconvolution and extrapolation to parton level 
with Monte Carlo simulations, and implementing the use of real observables as b-hadrons and 
b-jets. 



2 B hadron production cross section using J/ip 

A measurement of the B-hadron (H b ) production cross section has been performed by CDF, 
using data corresponding to an integrated luminosity of 39.7 pb _1 . First, the inclusive produc- 
tion cross section of J/tp is measured, considering J/i/j — > fi + pT . This cross section contains 
contributions from various sources, including prompt production of charmonium, decays of ex- 
cited charmonium states and decays of B-hadrons; once the latter J ftp component is extracted, 



it is possible to determine in a very clean way the B-hadron cross section. 

A di-muon trigger is used to select the events: two oppositely charged muons are required to 
be reconstructed from tracks measured in the tracking system and matched to track segments 
detected in the muon chambers. Central muons (|y M | <0.6) with p T >1.5 GeV/c are used: the 
offline reconstruction efficiency is about 99% indipendently by the pj,, thus the measurement of 
J/ip cross section covers all transverse momenta from to 20 GeV/c in the central rapidity range 
\y J /^\ <0.6. To separate prompt J/ip and J/ip from decay H b — > J/ipX, lifetime distributions 
are used. B-hadrons have long lifetimes, on the order of picoseconds, so that i/ip from those 
decays are usually displaced from the primary vertex, taken as the beam position in the r — <f> 
plane. The signed projection of the flight distance of J/tp on its transverse momentum, L xy , is, 
in general, a good measurement of the displaced vertex; nevertheless, to reduce the dependence 
on the J/ip p T bin, the variable x = L xy x (M J ^ /p^), called pseudo-proper decay time, is 
used. A Monte Carlo simulation is implemented to model the x distributions from B-hadron 
events and, finally, an unbinned maximum likelihood fit is used to extract the B fraction from 
the data for different ranges in Pj. . Figure ^ shows an example of fits to the pseudo-proper 
decay time; the B fraction varies from 10% (low p T ) to 45% (high p T ). The extracted differen- 
tial cross section of the B-hadrons over the transverse momentum range from to 25 GeV/c is 
shown in figure [21 The cross section is corrected for the branching fraction Br(H b — > J/tpX) 
= 1.16±0.10% and Br(J/i/j -> = 5.88±0.10%, and divided by two to have the single B- 

hadron cross section. The most recent theoretical calculations are also superimposed, showing 
a remarkable agreement with data. 
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Figure 1: Example of fits to the J/t/> pseudo- proper 
decay time in the range 1.25< p^F <1.5 GeV/c to 
extract J ftp from long-lived B-hadron decays. 
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Figure 2: B-hadron production cross section as a 
function of B-hadron p T . The crosses with erros 
bars are the data with systematic and statistical er- 
ror added. Theoretical prediction is superimposed. 



3 B-jet production 

The measurement of the inclusive 6-jet cross section^ in Runll (CDF), is based on about 300 
pb _1 of data. The use of jets extends the upper reach of the beauty production measurement 
performed using exclusive decays (p^ hadr < 25 GeV/c). Besides, theoretical uncertainties on 
fragmentation are rather small. A cone-based iterative MidPoint algorithm is used for jet recon- 
struction in the Y-cf) space, using a cone radius of 0.7. Only jets in the central rapidity region 
(|y| <0.7) are considered. The jet energy scale is corrected, using a Monte Carlo simulation, to 
compansate for energy losses at calorimeter level. 

The analysis exploits the good tracking capabilities of the detector and rely on b-jet identifi- 
cation done by secondary vertex reconstruction. The b-tagging algorithm uses displaced tracks 
associated with a jet that are within the jet cone. The search for secondary vertices is defined 
in two steps, with selection based on the significance of the impact parameter and of the decay 
lenght L . The sign of L is used as further selection to reject mistagged jets. To determine 




Figure 3: Left: Example of mass of secondary vertex distributions for data and fit prediction superimposed to 
Monte Carlo predictions, for one bin in jet P T (82< p J ^ 1 <90 GeV/c). Right: Fraction of b-tagged jets as a 
function of corrected jet P T : total errors and systematic uncertainties are superimposed. 

the heavy flavour content of a tagged jet, thus to extract the fraction of b-jets, the shape of the 
secondary vertex mass distribution is used: altough a full reconstruction of the hadron invariant 
mass is not possible, due to the presence of neutral particles and energy lost because of detector 
resolution, the invariant mass of tracks used to find the secondary vertex constitutes a good 
discrimination for jets from long-lived b or c hadrons. The mass of secondary vertex depends 
on jet transverse momentum and the fit is performed considering independently each bin in jet 
P T : in figure 01 on the left, mass distribution and fit results are presented for one bin in jet P T 
as example. On the right, the fraction of b-tagged jets as a function of jet P T is showen: the 
total errors and systematic uncertainties are superimposed. Figure |I] shows the inclusive b-jet 
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Figure 4: Left: b-jet cross section as a function of corrected jet P T , superimposed to Pythia Monte Carlo 
LO(CTEQ5L) prediction. Right: Ratio Data/Pythia as a function of corrected jet P T . 

cross section over a ~P T range between 38 and 400 GeV/c. Statistical and systematic error on 
the last bins is dominated by the error on b-tagged jets fraction. One of the main systematic 
error contributions comes from the jet energy scale: a conservative 5% uncertainty is considered 
at this stage, although on-going studies show the possibility to reduce it at 3%. The b-jet cross 
section is compared with leading order (LO) prevision from Pythia Monte Carlo (CTEQ5L), 
and results are in agreement with expectation: a NLO and beyond comparison is in progress. 
D0 is also performing a b-jet cross section measurement using a muon-tagging method to select 
jet sample enhanced in heavy flavor content: work is in progress to separate the b-jet content 
from the background. 



4 Upsilon production 

The D0 experiment has performed a precise measurement of the production of the T(1S) 
bottomonium stated, reconstructed through its decay Y(1S)— ► The differential cross 

section as a function of the T transverse momentum is determined in three rapidity ranges, 




Figure 5: Examples of di-muon mass distributions 
in rapidity range \y r \ <0.6 and 1.2< \y r \ <1.8, 
for range [4,6] GeV/c. The heavy line is the 
combined fit for signal and background. 



Figure 6: Top: Normalized differential cross sec- 
tions for T(1S) production. Theoretical prediction 
are superimposed. Bottom: Ratio of cross sections 
forcr(1.2< \y T \ <1.8) to cr(|y T | <0.6). Monte Carlo 
prediction is superimposed. 



< \y | < 0.6, 0.6< \y \ <1.2, and 1.2< \y \ <1.8. The data used in this analysis correspond 
to 159 pb _1 , collected with a scintillator-based di-muon trigger; two opposite-charge muon 
candidates are required, each of them with p T > 3 GeV/c and <2.2, and matched to a 
track found in the central tracking system. Isolation requirements are imposed to reduce the 
dominant background from bb. Two typical examples of di-muon mass distribution in different 
bins of rapidity and in the p T bin of 4 GeV/c< pj, <6 GeV/c are showed in [5J the clear 
T(1S) signal is accompained by a shoulder from T(2S) and T(3S). The mass distribution fits 
use three separate mass functions for each of the T(nS) states. The overall systematic errors 
are about 10%, excluding luminosity uncertainty that is 6.5%: the fitting procedure is one of 
the main source of uncertainties, while others as momentum resolution, track quality and track 
matching requirements contribute for less than 2% each. The current analysis assumes the 
T(1S) produced unpolarized, in agreement with previous CDF results (polarization parameter 
a=-0.12±0.22 for 8< p^ <20 GeV/c); thus, contribution to the systematic uncertainties is not 
included. Differential cross sections, normalized to unity, for the three rapidity ranges, are shown 
in figure El they are compared to theoretical predictions, that combined separate perturbative 
approacches for the low and high p T regions. 



5 Conclusions 



Studies of beauty production are important tests for perturbative QCD: results from CDF and 
D0 benefit from the very high statistic provided by the Tevatron collider and they are generally 
compatible with the recent theoretical predictions. Further improvements are expected in the 
near future and new measurements on bb cross section and correlations are already in progress. 



References 

1. S.Frixione, M.Mangano, N ucl. Phys. B 4 83 (1997) 321, |hep-ph/9605270| 

2. The CDF Collaborati on, |hep-ex/0412071[ Published in Phys.Rev.D71:032001 ,2005. 

3. CDF Public results in http://www-cdf.fnal.gov/physics/new/qcd/qcd_plots/ 
bjets_run2 /blessed_results.html 

4. TheD0 Collaboration, |hep-ex/ 0502030 , FERMILAB-PUB-05/020-E; Submitted to PRL. 



